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Transition State Theory forms the basis of computing reaction rates in chemical and other systems. 
Recently it has been shown how transition state theory can rigorously be realized in phase space using 
an explicit algorithm. The quantization has been demonstrated to lead to an efflcient procedure 
to compute cumulative reaction probabilities and the associated Gamov-Siegert resonances. In 
this letter these results are used to express the cumulative reaction probability as an absolutely 
convergent sum over periodic orbits contained in the transition state. 
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Introduction. — Transition State Theory, developed by 
Eyring, Polanyi and Wigner in the 1930's, is the most 
fundamental and widely used method to compute reac- 
tion rates. During a reaction a molecular system is en- 
visaged to pass through a 'transition state' or 'activated 
complex', a kind of unstable supermolecule poised be- 
tvifeen reactants and products [l[ . The main idea of tran- 
sition state theory is to place a dividing surface in the 
transition state region and compute the classical reac- 
tion rate from the directional flux through the dividing 
surface. In order not to overestimate the reaction rate 
the dividing surface needs to have the crucial property 
that it is crossed exactly once by all reactive trajectories 
(trajectories passing from reactants to products or vice 
versa) and not crossed at all by all other (non-reactive) 
trajectories. In the 1970's Pechukas, Pollak and others 
showed how to rigorously construct such a dividing sur- 
face from a periodic orbit giving the so-called periodic 
orbit dividing surface (PODS) [2|. The generalization to 
more degrees of freedoms has posed a major problem, and 
was solved only recently using ideas from dynamical sys- 
tems theory (see 0] ) . This shows that the transition state 
at energy E is formed by a normally hyperbolic invariant 
manifold (NHIM) (see Q), which in this case is an in- 
variant sphere of dimension 2(i— 3, where d is the number 
of degrees of freedoms, and normal hyperbolicity means 
that the contraction and expansion rates associated with 
the directions normal to the sphere dominate those of 
the directions tangential to the sphere. For d = 2, this 
simply is the unstable periodic orbit of the PODS [5|. In 
fact, the NHIM spans another sphere which is of dimen- 
sion 2d — 2 and hence has one dimension less than the 
energy surface and can be taken as a dividing surface. 
The NHIM forms the equator of this sphere and divides 
it into one hemisphere crossed exactly once by all forward 
reactive trajectories and one hemisphere crossed exactly 
once by all backward reactive trajectories. The NHIM 
itself is invariant and can be viewed as the energy sur- 
face of an invariant subsystem (the 'transition state' or 
'activated complex') with one degree of freedom less than 



the full system (i.e. with the reaction coordinate being 
frozen at a particular value) . All these phase space struc- 
tures can be explicitly constructed from a normal form 
which at the same time gives a simple expression for the 
flux through the dividing surface. In Q the quantization 
of this normal form has been used to develop a quantum 
version of transition state theory. This quantum normal 
form has been demonstrated to give an efficient method 
to compute cumulative reaction probabilities (the quan- 
tum analogue of the classical flux) and Gamov-Siegert 
resonances associated with the activated complex [a0]- 
In this letter we use these results to show that the cu- 
mulative reaction probability can be expressed as a sum 
over periodic orbits contained in the activated complex. 

The normal form representation of the activated com- 
plex and the computation of reaction rates. — We consider 
a molecular system with d = 1 + f degrees of freedom 
which has a saddle-center-. .. -center equilibrium point 
('saddle' for short), i.e. the matrix associated with the 
linearized Hamilton's equations has one pair of real eigen- 
values ±A, and / pairs of purely imaginary eigenvalues 
zbiwfc, k — 1,...,/. We will restrict ourselves to the 
generic case of linear frequencies ujk fulflUing no reso- 
nance condition miui + . . . -\-mfUjf — for any vector of 
integers m — (mi, . . . ,mf) ^ 0. Such saddles are char- 
acteristic for reaction type dynamics as for energies near 
the energy of the saddle, they induce a bottleneck type 
structure of the energy surface near the saddle through 
which the system has to pass in order to react. 

Normal form theory shows that in the neighbor- 
hood of the saddle there is a canonical transformation 
such that the transformed Hamiltonian is of the form 
IIo{I, Ji, . . . , Jf), where / = {Po ~ 9o)/2 is an inte- 
gral associated with the reaction coordinate, and the 
Jfc = {pI + ql)/2, k — 1,...,/, are action integrals 
associated with the bath modes. The activated com- 
plex is the invariant subsystem given by = =0. 
Its motions are described by the reduced Hamiltonian 
Hq{0, Ji, . . . , Jf), and thus is integrable, i.e. in action 
angle variables (J, ip) the equations of motion are J = 



2 



and ip = Vj-ffo(0, J) with solutions J(i) = const and 

ip{t) = ipa+tn{3) mod 2n , where J7(J) := VjiJo(0, J) . 

(1) 

The motion is thus quasiperiodic. It takesplace on invari- 
ant / dimensional Liouville- Arnold tori [8| which foliate 
the phase space of the activated complex. The motion 
becomes periodic for the J for which J1(J) — am, where 
m e Z-'^ and a e K,. We call the torus corresponding to 
this J a resonant torus. Fixing the energy E the energy 
surface of the activated complex, 



= {J G : i/o(0, Ji, . . . , J/) - E] , 



(2) 



is the action space projection of the NHIM mentioned 
in the introduction. The volume it encloses in the space 
of the actions J is proportional to the directional flux 
through the dividing surface (see Fig.[T]). 

A quantum normal form procedure based on the Weyl 
symbol calculus |6| shows that in the quantum mechan- 
ical case a unitary transformation can be found which 
transforms the Hamilton operator to the form H = 
H(I, Ji, . . . , Jn) which is a polynomial function of the 
operators / = {-h^dl^-ql)/2 and Jk = {-h'^d'^^ + ql)/2 
associated with the classical integrals. The polynomial 
defining the quantum normal form operator has the h 
expansion H{I,J) — Ho{I,3) + tiHi{I,3) -I- . . ., where 
Hk{I,J) are independent of H, and Ho{I,3) coincides 
with the classical normal form Hamiltonian. 

The cumulative reaction probability at energy E is 
then given by 
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where /„ — In{E) is implicitly defined by 

i?(/„,Ji,...,j/) = £;, 



(3) 



(4) 



and n = (ni, . . . ,nf) e INg is the vector of quantum 
numbers for the Bohr-Sommerfeld quantized actions, 

Jk = h{nk + Y)^ = (5) 

Here the ak = 2 are Maslov indices which for later refer- 
ence we group in the vector a — (ai, . . . ,af) (see Q for 
an earlier reference and @ where this result is derived in 
a systematic semiclassical expansion in h) . In the follow- 
ing we derive a formula which expresses N{E) in terms 
of a sum over periodic orbits. 

A periodic orbit formula for the cumulative reaction 
probability. — To derive our periodic orbit formula it is 
convenient to consider the energy derivative of the cumu- 
lative reaction probability (jS)), 
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FIG. 1: For d = 3 degrees of freedom, the left panel shows 
an energy surface Ho{I, J2, J3) ~ E for an energy above the 
saddle energy. The red lines mark the Bohr-Sommerfeld quan- 
tized actions J. The right panel shows the energy surface of 
the activated complex S_b defined in ([2| marked as the blue 
line in the left panel. The enclosed area is proportional to the 
classical fiux, and equivalently, to the mean number of states 
of the activated complex. 



Using (HI) the factor d/„/di? can be written as 



dE 



dH 



(7) 



We can obtain a periodic orbit formula for n{E) follow- 
ing a computation similar to the derivation of the Berry- 
Tabor trace formula for the density of states of classically 
integrable systems ■ Following we use the Poisson 
summation formula to rewrite ^ as 



27r 



n/2. 



dfj 



dH 
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i=i(E,j) 4cosh'{TTl{E,J)/n) 



where I{E, J) is determined by 

H{IiE,J),J,,...,Jf)^E. 



(8) 



(9) 



Note that the h expansion of the quantum normal form 
Hamiltonian implies, via (|9]), an h expansion of I{E,3), 
i.e. I{E, J) = Iq{E, J) + hh{E, J) 4^ . . .. In the following 
we separately discuss the term which we refer to as 
the Thomas- Fermi term [loj , and the remaining sum over 
m 7^ which we refer to as the oscillatory term nosc{E). 
The Thomas-Fermi term. — For m = 0, we get 
-1 



noiE) 



27r 



d^J 



dH 



This term can easily be interpreted from considering its 
integrated version 



No{E) 



dE'no(E') = ^ [ df J }^,^ . 

(11) 
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In the semiclassical limit, h 0, the integrand can be 
viewed as a characteristic function on the action space 
region I{E,3) > 0. The integral in ([TT|) hence gives the 
action space volume enclosed by the surface I{E, J) = 0, 
and accordingly Nq{E) is given by the classical flux di- 
vided by the elementary volume {2TTh)-^ , which agrees 
with the mean number of states of the activated complex 
to energy E @ (see Fig. [1]). The term riQ^E) is the cor- 
responding differential version, i.e. the mean density of 
states of the activated complex at energy E. 

The oscillatory term. — To compute the terms rim {E) 
for m 7^ we use 
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4cosh2(7rx) (27r)2 



, smh(?//2) 



(12) 



to rewrite ([5]) as 



n^{E) 



g— iTrmot /2 

27r;i/+i 



dy / d-^j( ^ I X 

J>0 V / I{E,3) 



y/2 



(13) 



i[2-Km3-yI(3,E)]/h 



sinh(2//2) 



This integral can be evaluated by the method of station- 
ary phase. The stationary phase conditions are 



2^m^yVjh{E,3) 



/o(i?,J)-0, (14) 



and by differentiating Hq{Iq{E^3),3) = E we obtain 



27rm = — 



A(J) 



n{J) , where A(J) 



dl 



(0,J), 



(15) 

The second condition in restricts J to the energy 
surface of the activated complex E e defined in ([2]) . The 
first conditions then fixes a point Jm on E^; (or a finite 
number of points Jm,i) by requiring that the frequency 
vector r2(J) at Jm is proportional to m. By ([T]) this 
means that the torus corresponding to Jm is resonant, 
and by we have 



(16) 



where y < (y > 0) if m and are parallel (anti- 
parallel). Here all functions of J are evaluated at Jm- 
Let Q be the (/-|- 1) x (/-I- 1) matrix of second derivatives 
of the phase function in evaluated at Jm and y, and 
/? its signature. We then find for n^i^E), 

(2^)^e-'["'"°/^+^"l'"l^^i/l"l+"^/4] y/2 



ft^A^ldetQI 



sinh(i//2) 



To evaluate the determinant of Q it is useful to intro- 
duce the curvature tensor K of E^. Let ei, 62, . . . , ey_i 
be / — 1 orthogonal unit vectors which are tangent to E^; 



at J. Noting that E^ is the hypersurface Io{E,J) 
we can write the components of K at Jm as 



1 



|Vj/o| 



ei ■ /g ej — - 



1 



iVjiJo 



(18) 



where /q and Hq denote the matrices of second deriva- 
tives with respect to J. Let ei be the unit vector parallel 
to Vji?o — Then in the basis of the Sj the matrix Q 
becomes 



\^i\/x 

\n\/X -yei/^ei 
a 



0^ 

y\VjIo\Kj 



(19) 



where a has components ei/gej. The determinant of 
this matrix can be evaluated straightforwardly, but to 
determine as well the signature it is useful to rewrite it 
as follows. Let A be the upper left 2x2 block of (fTOj) . 
D = y\VIo\K and B = (O a), then if det/C 7^ we can 
form 



A 

B D 



I B'^D-^ 
/ 



A- 



B^D-^B 
D 



I o^ 

D-^B I, 



By the special structure of B we find that B^D ^B = 
|j for some number c. Hence det(^ — B^D^^B) = 

-irjp/A^ < and so A- B'^D-^B has signature 0. The 
signature /? of Q is thus determined by D = y\WjIo\K 
and we find 

f3 — sign ?; sign if , (20) 
and with t/Vj/q = 27rm, by ((HI), the determinant is 



VidetQI = (27r|m|)^ V|detX| \n\/X 

evaluated at J = Jm- 

We notice that if Jm and y are a solution to the station- 
ary phase condition for m, then Jm and qy are a solution 
for qm for any q € Z\{0}. It is natural to choose /x ~ m 
with positive coprime components and combine the two 
terms with g/x and —q^- This way the n{E) contribution 
of the gth repetition of a resonant torus with ^ /x is 
given by 



2-K 



A 



ft(/+^)/^sinh(vr,MA) 
cos ((7(27r/2 - J/n - 7r/xQ;/2 - 27r|/x|A/i/|0|) + 7r;3/4) 



g(/-3)/2|/x|(/-3)/2|n|V|det/i(J'^)| 



(21) 



Example. — We consider the example of a system com- 
posed of an Eckart barrier and two Morse oscillators. Its 
quantum normal form Hamiltonian is given by ^] 




fc=i 



J 2m 



■Jk 



(22) 
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FIG. 2: Energy surfaces with resonance lines ^ii/ ^2 of the ac- 
tivated complex which consists of two Morse oscillators. The 
insets show the resonant tori with ^11,^2 < 3 projected to the 
configuration space of the oscillators. 



Here H — Hq, and hence Ii — 0. The frequencies are 

= dj.Ho = akJ— - ^ Jfe , = 1, 2 , (23) 

V TO TO 



and 



A 



dHo 



dl 



1=0 




(24) 



We choose D2 = ai — a2 — 1, -Di = 5/6, ag = 47r, 
Vb = 5/4, and h — 0.1. Figure [5] shows energy surfaces 
of the activated complex consisting of the two Morse 
oscillators together with some resonance lines rii/f22 = 
Ah/M2. The sign of the curvature matrix is = — 1. The 
exact cumulative reaction probability, and its derivative 
n{E) can be computed analytically for this system @. Its 
oscillatory part, tIosc = ^[E) — no{E), is shown together 
with its approximation by the periodic orbit sum over 
the terms (f2T|) for jUi , /X2 < 3 in Fig. [H 

Conclusions. — In this letter we derived a periodic or- 
bit formula for the cumulative reaction probability, and 
demonstrated its applicability for a simple example. In 
the limit A — (no tunneling through the potential bar- 
rier) our periodic orbit formula reduces to the Berry- 
Tabor trace formula for the density of states of the acti- 
vated complex. In the general case A 7^ 0, our periodic 
orbit formula is (as opposed to the Berry- Tabor trace for- 
mula) absolutely convergent due to an additional factor 
which leads to an exponential damping of contributions 
of long periodic orbits. Although we incorporated only 
six periodic obits (and their repetitions) in our example 
the agreement with the exact result is already very good. 
This is even more impressive as we have so far only taken 




FIG. 3: Exact (dashed line) and periodic orbit approximation 
(solid line) of the energy derivative of the cumulative reaction 
probability including resonant tori with /ii , /i2 < 3. 



into account simple stationary points associated with res- 
onant tori, and no isolated and ghost orbits which would 
naturally arise in a more elaborate uniform approxima- 
tion [13, • Similarly, the integral associated with the 
reaction direction can be cast into a periodic orbit sum 
over the instanton orbits Q extending the applicability 
of our periodic orbit formula to energies below the saddle 
energy. These aspects will be discussed in more detail in 
a longer version of this letter. 
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